ORIGINAL ARTICLE 

http://dx.doi.org/1 0.4306/pi.201 3.1 0.2.1 73 



Print ISSN 1738-3684 / On-line ISSN 1976-3026 

OPEN ACCESS 



Impact of the BDNF Val66Met Polymorphism on Regional Brain 
Gray Matter Volumes: Relevance to the Stress Response 



Sung Nyun Kim 1,3 , Do-Hyung Kang 1,4 , Je-Yeon Yun 14 , Tae Young Lee 4 , 
Wi Hoon Jung 3 , Joon Hwan Jang 1,4 and Jun Soo Kwon 1,2 ' 3,43 

'Department of Psychiatry, Seoul National University College of Medicine, Seoul, Republic of Korea 

department of Brain and Cognitive Sciences - World Class University Program, Seoul National University College of Natural Sciences, Seoul, 
Republic of Korea 

institute of Human Behavioral Medicine, SNU-MRC, Seoul, Republic of Korea 
"Department of Neuropsychiatry, Seoul National University Hospital, Seoul, Republic of Korea 



Objective Genetic imaging is used to investigate the mechanism by which genetic variants influence brain structure. In a previous 
study, a structural change of the dorsolateral prefrontal cortex was associated with symptom modulation in post-traumatic stress disor- 
der patients. This study examined the effect of a polymorphism in the gene encoding brain-derived neurotrophic factor (BDNF) on re- 
gional gray matter (GM) volumes and the correlations between the dorsolateral prefrontal GM volume and the stress level in healthy 
volunteers. 

Methods Sixty-one volunteers underwent genotyping for the BDNF Val66Met single nucleotide polymorphism (SNP) and completed 
the Stress Response Inventory (SRI). Magnetic resonance images were also acquired, and the effect of each subjects BDNF genotype and 
SRI subscore on his or her dorsolateral prefrontal GM volume was evaluated. 

Results The Val/Val homozygotes had significantly larger GM volumes in the prefrontal cortex and the precuneus, the uncus, and the 
superior temporal and occipital cortices than Met carriers. The Met homozygotes demonstrated a higher stress response in depression 
domain than Val/Val and Val/Met groups. A negative correlation between the middle frontal cortex GM volume and the SRI depression 
subscore was found. 

Conclusion These findings indicate an interaction between genes and brain structure, and they suggest that differences in dorsolateral 
prefrontal GM volume related to the BDNF Val66Met SNP are associated with resilience to stressful life events, particularly in the di- 
mension of emotion. Psychiatry Investig 2013;10:173-179 

Key Words Brain-derived neurotrophic factor, Val66Met, Genetic polymorphism, Magnetic resonance imaging, Stress, 
Dorsolateral prefrontal cortex. 



INTRODUCTION ticity in response to stressful life events. Acute and chronic 
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Brain-derived neurotrophic factor (BDNF), a neurotroph- 
in, plays essential roles in neuronal proliferation, differentia- 
tion, and survival during development and functions in the 
molecular mechanisms underlying neural plasticity and con- 
nectivity. In addition, BDNF is associated with neuronal plas- 



stresses increase BDNF expression in the basolateral amygda- 
la. 1 Moreover, BDNF itself stimulates hypothalamic hormone 
synthesis, 2 which might be directly involved in the physiologi- 
cal response to stress. Additionally, BDNF-related stress resil- 
ience is conducted via the mesolimbic pathways, which are 
associated with bidirectional projections between the ventral 
tegmental area versus the nucleus accumbens and the pre- 
frontal cortex (PFC). 3 The PFC is a region of abundant expre- 
ssion among the many areas of the central nervous system in 
which BDNF is expressed. 



An important polymorphism has been identified in codon 
66 of the gene encoding BDNF, which is located on chromo- 
some llpl4.1. This single nucleotide polymorphism (SNP) 
causes the substitution of methionine (Met) for valine (Val) 
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(Val66Met). The Met allele has been linked with functional 
impairment of the hippocampus, including that of working 
memory performance. 4 Additionally, the presence of the Met 
allele is associated with a decreased gray matter (GM) volume 
throughout the PFC, including the anterior cingulate cortex 
(ACC), 5 ' 6 the dorsolateral PFC (DLPFC), 7 and the medial tem- 
poral lobes, which include the hippocampus, parahippocam- 
pal gyrus, and amygdala. 8 

Several recent studies have suggested that gene-stress inter- 
actions are mediated by regional brain volume differences ac- 
cording to the SNP allelic composition; the BDNF Val66Met 
SNP and DLPFC may be relevant candidates. In previous 
studies of subjects with post-traumatic stress disorder (PTSD) 
who had experienced traumatic stress, the structure and func- 
tion of the DLPFC were reported to be altered; these changes 
were hypothesized to play an important role in modulating 
the fear and anxiety circuits. 9 A recent neuroimaging study 
showed that greater DLPFC thickness is associated with 
greater PTSD symptom improvement and that the BDNF 
Val66Met gene polymorphism is associated with DLPFC thick- 
ness in PTSD patients. 10 The authors concluded that the DLP- 
FC thickness increased via the experience-dependent synap- 
tic plasticity mediated by BDNF genetic variation in PTSD 
patients. However, the roles BDNF and the DLPFC play in 
mediating stress in healthy subjects remain unclear. 

The present study examined regional GM volume differ- 
ences in relation to the BDNF Val66Met SNP in normal sub- 
jects. Additionally, regarding the hypothesized role of the 
DLPFC in corticosteroid-related stress reactivity, which is un- 
der BDNF-related physiologic control, we analyzed the corre- 
lation between the DLPFC GM volume and the level of stress 
in 61 healthy volunteers. We sought to elucidate the structural 
neural correlate of the stress response in relation to the BDNF 
Val66Met SNP in the normal population. 

METHODS 
Participants 

A total of 61 healthy volunteers (35 males and 26 females) 
were recruited as subjects and were evaluated using the Struc- 
tured Clinical Interview for DSM-IV non-patient version. 
None of the subjects had been diagnosed with psychosis, bi- 
polar disorder, major depressive disorder, substance abuse or 
dependence, significant head injury, seizure disorder, or an 
IQ lower than 70. This study was approved by the Institution- 
al Review Board at Seoul National University Hospital, and 
written informed consent was obtained from all participants 
after the procedures had been fully explained. This study was 
conducted in accordance with the Declaration of Helsinki. 

All participants completed the Stress Response Inventory 



(SRI), 11 a self-report questionnaire that evaluates stress re- 
sponses. The SRI includes 22 questions designed to address 
the following dimensions: somatization, depression, and an- 
ger. Each question was scored on a 5-point Likert scale, from 
0 ("not at all") to 4 ("absolutely"). The reliability of the SRI 
was confirmed by measuring Cronbachs alpha (somatization, 
0.89; depression, 0.88; and anger, 0.87). 12 

Genotyping of BDNF 

The subjects' DNA was genotyped for the BDNF Val66Met 
SNP. Whole blood (12 mL) was drawn from each participant 
by venipuncture into EDTA-containing tubes. After measur- 
ing the blood volume, a 1-mL aliquot was stored at -20°C un- 
til genomic DNA extraction. Genomic DNA was isolated from 
200 uL of whole blood using a QIAamp DNA Blood Mini Kit 
(Qiagen, Hilden, Germany). The BDNF polymorphism was 
detected by TaqMan allelic discrimination assays using an 
ABI Prism 7500 Sequence Detection System (Applied Biosys- 
tems, Foster City, CA, USA), by which the BDNF Val66Met 
(rs6265, C_11592758_10) SNP was investigated. The reac- 
tion mixture contained 5 uL of 2x TaqMan genotyping mas- 
ter mix, 0.5 uL of 20 x drug metabolism genotyping assay 
mix, 3.5 uL of DNase-free water, and 1 uL of genomic DNA 
in a total volume of 10 uL. The products were analyzed using 
ABI Prism 7500 Sequence Detection System software (version 
2.0.1). 

Magnetic resonance image acquisition and 
preprocessing 

Magnetic resonance images were acquired on a 1.5 T Mag- 
netom Avanto scanner (Siemens, Erlangen, Germany) using a 
Tl-weighted magnetization -prepared rapid acquisition gradi- 
ent echo sequence with the following imaging parameters: ac- 
quisition matrix=256x256, voxel size=0.45x0.45x0.90 mm 3 , 
axial slices=176, field of view=23x23 cm 2 , repetition time= 
1,160 ms, echo time=4.76 ms, and no slice gap. 

The raw DICOM magnetic resonance imaging (MRI) data 
were converted to the NIFTI format using SPM8 software 
(http://www.fil.ion.ucl.ac.uk/spm/software/spm8/) running 
on MATLAB R2010B (The Math Works Inc., Natick, MA, 
USA). Brain volume normalizing, bias correcting, and seg- 
mentation into GM, white matter (WM), and cerebrospinal 
fluid (CSF) were performed with the VBM8 toolbox (http:// 
dbm.neuro.uni-jena.de/vbm/), which is incorporated in 
SPM8. This method uses an optimized voxel-based mor- 
phometry (VBM) protocol and a model based on Hidden 
Markov Random Fields, which was developed to increase the 
signal-to-noise ratio. The total volumes of GM, WM, and CSF 
were calculated by adding the resulting tissue probabilities. 
Total brain volume was defined as the sum of the GM and 
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Table 1 . Demographic data, genetic polymorphism profile of subjects and the SRI scores 





Total 


BDNF Val/Val 


BDNF Val/Met 


BDNF Met/Met 








(N=61) 


(N=20) 


(N=29) 


(N=12) 


F 


P 


Ase 


25.79±4.01 


24.95±4.55 


26.52±3.75 


25.42±3.66 


0.967 


0.386 


Sex (male/female) 


35/26 


14/6 


16/13 


5/7 


2.530 


0.282 


IQ 


114.07±10.26 


116.79±10.72 


113.82±7.96 


110.33±13.53 


1.497 


0.233 


SRI total score 


16.84±12.54 


13.85±10.69 


16.97± 13.00 


21.50±13.77 


1.418 


0.250 


SRI somatization 


6.6615.16 


5.70±4.78 


7.31 ±5.50 


6.67±5.09 


0.568 


0.570 


SRI depression 


6.49±6.17 


4.5515.20 


6.00±5.56 


10.91±7.30 


4.682 


0.013* 


SRI anger 


2.13±1.94 


2. 11 ±2.02 


2.00±1.79 


2.50±2.24 


0.280 


0.757 



Means±SD. *p<0.05/3=0.017. BDNF: brain-derived neurotrophic factor, Val: valine, Met: methionine, SRI: stress response inventory 



WM volumes. Afterwards, the volume images were smoothed 
using an isotropic Gaussian kernel (full width at half maxi- 
mum=8 mm). 

Statistical analysis 

Chi-square tests were performed to verify Hardy- Weinberg 
equilibrium and to compare the allelic and genotypic frequ- 
encies. 

Analysis of variance (ANOVA) was used to compare de- 
mographic data and three SRI subscores (somatization, de- 
pression, and anger) among the Val/Val, Val/Met, and Met/ 
Met groups, and a post-hoc analysis was conducted for pair- 
wise comparisons between each genotype group. 

Voxelwise comparisons was applied to the GM images us- 
ing genotype (Val/Val homozygotes vs Met+ carriers) as the 
factor and age and sex as covariates. All statistical tests for 
GM imaging were executed with a threshold of uncorrected 
p<0.005 and k>50. For the region of interest (ROI) analysis, 
statistically significant clusters showing the greater brain vol- 
ume of the Val/Val group compared to the Met+ group were 
extracted using the MarsBaR toolbox (version 0.42). 

Correlations between the SRI subscores and ROI cluster of 
the dorsolateral prefrontal cortex were also investigated. The 
significance threshold was set at p<0.05, and a Bonferroni 
correction (p<0.017) was applied to the SRI subscores ANO- 
VA and correlation analysis. Tukey's test was used for a post- 
hoc pairwise analysis. All analyses were conducted using 
PASW version 18.0 (Statistical Package for the Social Scienc- 
es, SPSS Inc., Chicago, IL, USA) and SPM8 (Statistical Para- 
metric Mapping, Wellcome Department, University College 
of London, UK, http://www.fil.ion.ucl.ac.uk/spm). 

RESULTS 

Demographic and genetic characteristics 

Our study included 61 subjects in total and grouped them 
according to BDNF genotype into 20 Val allele homozygotes 




Figure 1. Significant differences in gray matter volume in relation 
to the BDNF SNP projected on a standard brain template (Signifi- 
cance thresholds: uncorrected p<0.005 and k>50). BDNF: brain- 
derived neurotrophic factor, Val: valine, Met: methionine, SNP: 
single nucleotide polymorphism. 



(Val/Val), 29 heterozygotes (Val/Met), and 12 Met allele ho- 
mozygotes (Met/Met). The frequency of the BDNF genotypes 
was 32.8% for Val/Val, 47.5% for Val/Met, and 19.7% for Met/ 
Met in our sample. There were 20 Val homozygotes and 41 
Met allele carriers. There were no statistically significant dif- 
ferences in age, sex, or IQ between the genotype groups (Table 
1). No significant deviation from Hardy- Weinberg equilibri- 
um was found (x 2 =0.06, p>0.05). 
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Val/Val Val/Met Met/Met 



Figure 2. Significant differences in SRI depression scores accord- 
ing to BDNF SNP genotypes. *p<0.05, multiple comparison cor- 
rected. SRI: stress response inventory, BDNF: brain-derived neu- 
rotrophic factor, SNP: single nucleotide polymorphism, Val: valine, 
Met: methionine. 

Effect of BDNF genotype on regional brain 
GM volumes 

VBM revealed regional brain GM differences between the 
two groups (Val/Val and Met+)(Table 2). In the DLPFC re- 
gion, left inferior frontal gyrus (BA 45) and middle frontal 
gyrus (BA 9) were smaller in Met-I- group. Met+ carriers also 
had smaller GM volumes in the right culmen of the cerebel- 
lum, right uncus of the basal ganglia (BA 36 and 28), right su- 
perior temporal gyrus (BA 38), left superior occipital gyrus 
(BA 19), left precuneus (BA 19), and left superior parietal 
lobule (BA 7) than Val homozygotes. In contrast, the regional 
GM volumes of the left superior frontal gyrus (BA 10) and 
left superior parietal lobule (BA 7) were larger in the Met+ 
group (uncorrected p<0.005, k>50)(Figure 1). 

Effect of BDNF genotype on the stress response 

Groups of different BDNF genotypes showed a significant 
difference on the SRI depression subscore (F=4.682, p=0.013) 
(Figure 2). Met/Met homozygotes had a higher SRI depres- 
sion subscore than Val/Val homozygotes (p=0.011, Tukey 
corrected) and Val/Met heterozygotes (p=0.044, Tukey cor- 
rected) in the post-hoc analysis. Regarding the somatization 
and anger SRI scores, no statistically significant differences 
were found according to the genotypes. 

Correlations between the SRI subscores and ROIs 

We selected a significant cluster from the left inferior fron- 
tal and middle frontal cortex as ROIs to further investigate 
the correlations between the three SRI subscores (somatiza- 
tion, depression, and anger) and GM volume of the DLPFC. 
A statistically significant correlation was found between the 
SRI depression subscore and the left middle frontal GM vol- 
ume for all 61 subjects [Pearson's correlation coefficient= 
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Figure 3. Scatterplot depicting the distributions of the SRI depres- 
sion subscore and the brain GM volume of the left middle frontal 
gyrus in all subjects. Solid line: trend of the correlation between the 
two variables. SRI: stress response inventory, GM: gray matter, 
Val: valine, Met: methionine. 



-0.313, p=0.014; p<0.05/3 (Bonferroni corrected)] (Figure 3, 
Table 3). There were no significant correlations between the 
inferior frontal gyrus GM volume and the SRI subscores. 

DISCUSSION 

In this study, we aimed to elucidate the structural neural 
correlate of stress resilience in relation to the BDNF Val66Met 
polymorphism in normal subjects. We found that Met/Met 
homozygotes had a higher stress response in the depression 
domain (SRI depression subscore) compared to the Val/Val 
and Val/Met groups. By investigating whole-brain GM vol- 
ume differences according to the BDNF Val66Met SNP, we 
found that brain areas such as the middle frontal gyrus, the 
inferior frontal gyrus, precuneus, uncus, and superior tempo- 
ral cortex had smaller GM volumes in the BDNF Met carriers 
compared to Val homozygotes. Additionally, we found a neg- 
ative correlation between the GM volume of the middle fron- 
tal cortex (BA 9) and the stress response measured by the SRI 
depression subscore. To our knowledge, this is one of the first 
reports to demonstrate that the BDNF Val66Met SNP medi- 
ates resilience to stressful life events by influencing regional 
GM volume, particularly in the DLPFC region. 

The middle and inferior frontal cortices exhibited greater 
GM volumes in the Val homozygotes compared to the Met 
allele carriers in this study. This result is in line with those of a 
previous VBM study targeting patients with bipolar disorder, 
which reported smaller GM volumes in the DLPFC and hip- 
pocampus in Val/Met bipolar patients compared with Val ho- 
mozygote patients. 13 In addition, the presence of the Met al- 
lele was reported to be associated with a more profound and 
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Table 2. Regional differences in brain gray matter volume between BDNF Val homozygote group (N=20) and BDNF Met+group (N=41)* 



MNI coordinates 

Brain area Laterality BA T Z Cluster size p value 

x y z 



Val/Val>Met+ 



Middle frontal gyrus 


L 


9 


-24 


15 


36 


3.29 


3.14 


72 


0.001 


Triangular part of inferior frontal gyrus 


L 


45 


-30 


15 


22 


3.35 


3.18 


124 


0.001 


Superior occipital gyrus 


L 


19 


-33 


-72 


22 


3.44 


3.27 


575 


0.001 


Precuneus 


L 


19 


-29 


-78 


39 


3.34 


3.18 




0.001 


Superior parietal lobule 


L 


7 


-41 


-75 


49 


3.23 


3.08 




0.001 


Oilmen of cerebellum 


R 




18 


-43 


-8 


3.81 


3.58 


250 


<0.001 


Uncus 


R 


36 


26 


2 


-33 


3.52 


3.34 


427 


<0.001 


Uncus 


R 


28 


17 


5 


-33 


3.23 


3.08 




0.001 


Superior temporal gyrus 


R 


38 


26 


15 


-33 


3.05 


2.93 




0.002 


Met+>Val/Val 




















Superior parietal lobule 


L 


7 


-29 


-46 


58 


3.29 


3.13 


88 


0.001 


Superior frontal gyrus (OFC) 


L 


10 


-24 


56 


3 


3.27 


3.12 


273 


0.001 


Superior frontal gyrus 


L 


10 


-20 


66 


-2 


3.22 


3.08 




0.001 



*p<0.005 and k>50. BDNF: brain-derived neurotrophic factor, Val: valine, Met: methionine, BA: Brodmann area, MNI: Montreal Neurologi- 
cal Institute, OFC: orbitofrontal cortex 



Table 3. Partial correlation analyses showing degrees of correlation between SRI and gray matter volume of left middle frontal gyrus in the 
whole subjects (N=61 ) 

SRI somatization SRI depression SRI anxiety 

GM volume of LtMFG Pearsons correlation coefficient -0 296 -0 313 -0 094 

p value (two-tailed) 0 020 0.014* 0.473 

*p<0.05/3=0.017. SRI: stress response inventory, GM: gray matter, MFG: middle frontal gyrus 



widespread reduction in the bilateral DLPFC GM volume in 
a healthy aged population. 14 The DLPFC has been identified 
as the neural correlate of psychological recovery from a se- 
verely stressful event because the degree of cortical thickness 
enrichment in the DLPFC area after trauma exposure could 
successfully predict the reduction in PTSD symptoms. 10 In 
particular, we showed a correlation between a greater middle 
frontal gyrus (BA 9) GM volume and the lower SRI depres- 
sion subscore. This result is in line with a previous functional 
MRI (fMRI) study of acutely traumatized individuals using a 
script-driven trauma imagery paradigm, in which BOLD ac- 
tivation of the middle frontal cortex (BA 47) was correlated 
with resilience to the stressful event. 15 In another fMRI study 
using the emotional face paradigm, DLPFC hypo-activation 
was significantly correlated with the degree of cognitive vul- 
nerability to stress, in-line with the diathesis-stress theory of 
depression. 16 

We found that several posterior brain regions were corre- 
lated with the BDNF SNP and could function as neural corre- 
lates of various stress reactions. In the 99m-HMPAO SPECT 
study, which compared patients with PTSD and trauma-ex- 
posed non-symptomatic controls, differences in HMPAO 



tracer distribution were found in the uncus. 17 Likewise, recall 
of trauma-related autobiographic memory in PTSD patients 
was related to increased cerebral blood flow in the precune- 
us, 18 which was negatively correlated with the WHO-10 score, 
a measure of current subjective well-being. In the temporal 
area, cortical thinning in the right superior temporal gyrus of 
subjects with PTSD was associated with an M50 auditory sen- 
sory gating deficit. 19 The superior temporal and occipital (BA 
19) cortices also showed increased activation in traumatic sc- 
ript-driven dissociative responses in PTSD patients who had 
been sexually abused. 20 

The superior parietal lobule (SPL) is regarded as a compo- 
nent of the mirror neuron system and also as a neural corre- 
late of analytical reasoning, including cognitive coping strate- 
gies for stressful life events. 21,22 In this study, a group compa- 
rison according to the BDNF Val66Met SNP revealed diverse 
GM thicknesses along the rostro-caudal axis of the left SPL 
(BA 7). Specifically, the GM volume of the posterior portion 
of the SPL (7P) was thicker in the Val homozygote group, wh- 
ereas that of the SPL anterior portion (7A) was thicker in the 
Met+ group. The left SPL area showed decreased activation 
during exposure to an emotional stress paradigm in a previ- 
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ous fMRI study of somatoform pain disorder patients. 23 This 
result needs further investigation. 

This study had several limitations that should be consid- 
ered. First, we did not consider the presence and intensity of 
previous stressful life events, despite the possible effect of 
BDNF expression. Likewise, physiologic responses to stressful 
events, including increased Cortisol levels, could have a delete- 
rious effect on regional brain volumes, including those of the 
hippocampus and prefrontal cortices. 24,25 Second, we had to 
merge the Val/Met and Met/Met groups as Met carriers in the 
GM analysis because of the small number of Met homozy- 
gotes (n=12). Recently, an MRI study of Met/Met individuals 
argued that they had a wider GM difference than Val/Met 
heterozygotes compared to Val/Val homozygotes. 26 Due to 
the small allele frequency of Met, further research is needed 
with a larger sample size. 

In this VBM-based study of regional GM volume differ- 
ences according to the BDNF Val66Met SNP, we found that 
the dorsolateral prefrontal region and multiple posterior brain 
regions in the cerebellum, visual association cortex, superior 
temporal cortex, uncus, and precuneus showed significantly 
decreased GM volumes in the Met+ carriers compared to Val 
homozygotes. Met/Met homozygotes demonstrated a higher 
stress response in the depression domain compared to the 
Val/Val and Val/Met groups. Moreover, a negative correlation 
was found between the GM volume of the DLPFC and the 
SRI depression subscore, suggesting that regional GM volume 
differences according to the BDNF Val66Met SNP are associ- 
ated with resilience to stressful life events, especially in the di- 
mension of mood regulation. A future study exploring the in- 
teraction between BDNF expression, stressful life events, 
regional brain volumes, and resilience to stressful events sh- 
ould help elucidate these complex gene-environment interac- 
tions. 
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